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Abstract 

We consider a standard susceptible-infected-susceptible (SIS) model to 
study behaviors of steady states of epidemic spreading in small-world net- 
works. Using analytical methods and large scale simulations, we recover the 
usual epidemic behavior with a critical threshold A c below which infectious 
diseases die out. Whereas for the spreading rate A far above A c , it was found 
that the density of infected individuals p as a function of A has the property 
p«F(iv")(lnA-lnA c ). 
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Complex networks have attracted an increasing interest recently. The main rea- 
son is that they play an important role in the understanding of complex behaviors 
in real world net works, JU El including the structure of language,El|3 scientific 
collaboration networks, the Internetd and World Wide Web,[0 E3 power 
grids, HIl food webs. lfl2l fT3ll chemical reaction networks . ifTH metabolic[15] and 
protein networks. lfT6l etc. In the study of complex networks, an analysis of the struc- 
tures can give important information about the underlying processes responsible for 
the observed macroscopic behavior, such as small-world lfTTl and scale-free[ 18 1 net- 
works. In particular, social networks have two characters. First, they show "cluster- 
ing", meaning that two of your friends are far more likely also to be friends of each 
other than two people chosen from the population at random. Second, they exhibit 
the "small-world effect", namely, that any two people can establish contact by go- 
ing through only a short chain of intermediate acquaintances. These two properties 
appear contradictory because the first is a typical property of low-dimensional lat- 
tices but not of random graphs or other high-dimensional lattices, while the second 
is typically of random graphs, but not of low-dimensional lattices. D.J. Watts and 
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S.H. Strogatz suggested a new small- world model[17| recently which interpolates 
between low-dimensional lattices and random graphs and displays the both proper- 
ties. In the model, a small- world network is constructed as follows: starting with a 
ring of N vertices, each connected to its 2K nearest neighbors by undirected edges, 
and then each local link is visited once with the rewiring probability p it is removed 
and reconnected to a randomly chosen node. Duplicate edges are forbidden. After 
the whole sweep of the entire network, a small-world network is constructed with 
an average connectivity (k) = 2K. The WS networks has been widely studied be- 
cause it constitutes an interesting attempt to translate complex topologies of social, 
economic, and physical networks into a simple model. 

Although topological properties of complex networks have been studied in detail. lfT8l 
[l9j|20l|2Tl a natural question arises, that is the dynamical properties which result 
from different networks. A good example is to inspect complex features of epidemic 
spreading since the characterization and understanding of epidemic dynamics in 
these networks could probably provide us immediate applications to a large number 
of problems, such as computer virus infections, distribution of wealth, transmission 
of public opinion, etc. Recent papers ll22l l23l l24l [231 have given some valuable in- 
sights of that: for small-world networks, there is a critical threshold below which 
an infection with a spreading rate dies out; on the contrary, for scale-free networks, 
even an infection with a low spreading rate will prevalence the entire population. 

In the present paper we consider a standard SIS model[26] in small- world net- 
works, in which each node represents an individual of the population and edges 
represent physical interactions through which an infection spreads. According to 
the SIS model, an individual is described by a single dynamical variable adopting 
one of the two stages: susceptible and infected. A susceptible individual at time 
t — 1 will pass to the infected state with the rate v at time t if it is connected to one 
or more infected individuals. Infected individuals at time t — 1 will pass to the sus- 
ceptible state with the rate 5 at time t, defining an effective spreading rate A = v/5. 
We can still keep generality by setting 5=1. Individuals run stochastically through 
the cycle susceptible — > infected — > susceptible, hence the model got its name. In 
the SIS model, an important observable is the prevalence p, which is the time aver- 
age of the fraction of infected individuals reached after a transient from the initial 
condition. Given a network, the only parameter of the model is the spreading rate A. 
The information on the global spreading of the infection is contained in the function 
p(A). 

Small-world networks have very small diameters which mean the presence of 
disordered long range interactions. In this case, the networks are very homogeneous 
and quite reasonable that the mean-field (MF) method is valid. By neglecting the 
density correlations among the different nodes and ignoring all higher order correc- 
tions in p(t), the time evolution equation of the SIS model can be written as ll2"7l 



2 



P(t) = -P(t) + Hk)p(t)(l-P(t)). (1) 

In the equation, the first term on the right-hand side (rhs) considers infected nodes 
become healthy with unit rate and the second term on the rhs represents the average 
density of newly infected nodes generated by susceptible nodes. By imposing the 
stationary condition d t p(t) = 0, one can obtain the equation 

p[-l + X(k)(l-p)] = Q (2) 

for the steady state density of infected nodes p. The equation defines an epidemic 
threshold 

A c HA;)- 1 . (3) 

In other words, if the value of A is above the threshold, A > A c , the infection spreads 
and becomes endemic with a finite stationary density p. Below it, A < A c , the 
infection dies out. In Euclidean lattices, J. Marro and R. Dickman have concluded 
the order parameter behavior of critical phenomena is p ~ (A — Ac)' 3 with (3 < 1 in 
the region A ~ A c ,|27| which uncovers the linear property in the critical dimension. 
R. Pastor-Satorras and A. Vespignani also recovered this property in small-world 
networks with the rewiring probability p = 1.0 [24J recently; it is worth noticing 
that in the extreme case the generated network is an entirely random network with a 
restriction which leads to a large cluster. 

In order to compare with the analytical prediction we have performed large scale 
simulations of the SIS model with parallel updating in the WS network with the 
rewiring probability p = 0.1. The both properties of social networks are well pre- 
sented by the network in this case. The size of the network is N = 10 6 . The number 
of the initially infected nodes is 10 percent of the size of the network. Simulations 
were implemented on the network averaging over 20 different realizations. After an 
initial transient regime, the systems stabilize in a steady state with a constant aver- 
age density of infected nodes. In Fig. [l]we plot the steady density of infected nodes 
p with various types of the scale of axis for A is very closed to A c . The linear prop- 
erty of the order parameter of Euclidean lattices, p ~ (A — Ac)' 3 , is well presented by 
the log-log plot of Fig. [l{d) which gives the parameter value (3 = 0.98 ± 0.04. How- 
ever, it is obvious that the four plots almost show the same shapes. Consequently 
the following acceptable predictions for Fig. l(a-d) are respectively given by 



p m ai\ + bi (4a) 

p « a 2 In A + b 2 (4b) 

p « a 3 (A - A c ) + h (4c) 

In p « a 4 ln(A — A c ) + b 4 (4d) 
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Figure 1: Density of infected nodes p in the WS network with K = 3 from the 
simulations near by A c = 0.1708 ± 0.004, which is in good agreement with the MF 
predictions A c = 1/2K = 0.1667. The numerical results were presented with p vs 
A (a), p vs log A (b), p vs A — A c (c), and logp vs log (A — A c ) (d). It indicates that 
all the plots perform the similar behaviors. 

Considering the Taylor expansion In A ~ j- + In A c — 1 at one order near A c , Eq. 
(PEEl) can be rewritten as 

p « a' 2 X + b' 2 . (5) 

One can also derive the equation 

p « <(A - A c ) + b' A (6) 

from the Eq. d4dh in the same way. So, based on the numerical simulations closed 
to A c , a analogous linear relationship between the infected density and the spreading 
rate, p ~ a\ + b, was extracted from the simple approximation of series expansions 
(see Eq. d4al) . Eq. ©, Eq. d4cT) . and Eq. © corresponding to Fig. da)-(d), 
respectively). Therefore, we will calculate the tendency of p for A > A c next. 

As explicitly stated in Ref. E71 . Eq. ^ is derived in the case of small values 
of p. So it is meaningful only in the region A ~ A c and could not determine the 
behavior for A far above A c , i.e., the MF prediction of Euclidean lattices is invalid 
in this region. In this condition, numerical methods are naturally adopted to find the 
behavior of the prevalence for A > A c . In Fig. |2]we perform the numerical results as 
far as A = 0.4. The linear property of Eq. ( PfaT ) was presented excellent agreement 
with the numerical results. Note that the constant parameters a 4 and 6 4 in Eq. ( PEdT ) 
should be taken count of the initial setting value of the average connectivity 2K. 
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Figure 2: Extensive numerical results of the density of infected nodes p in the WS 
network for A > A c . All plots perform the identical behavior which is described by 
Eq. ©. Parameter values (from right to left) K — 3, 4, 5. 

More exact conclusion for the nature of Fig. El which presented the whole behavior 
at A > A c , should be given by 

p^F(K) (In A - In A c ) , if A > A c . (7) 

To complete our study of the steady states of the SIS model in small- world networks, 
we compute the coefficient F(K) in Eq. ©. In Fig. |3]we plot the coefficient F(K) 
as a function of K. The linear behavior of the log-log plot predicts the power law, 
F = K~ a . Numerical results give the parameter value a ~ 0.98. 

In summary, we have analytically and numerically studied the steady states of 
epidemic spreading in small-world networks. In the region A ~ A c , we recover the 
MF results of Euclidean lattices, p ~ (A — A c ) s . But when spreading rates become 
far above the threshold, A > A c , the MF method can not work normally and the 
numerical method is adopted. It was found that the behavior of order parameter has 
the property p m F(K)(\n\ — In A c ). In the present work, we just consider the 
dynamics on the WS network with the rewiring probability p = 0.1. However, for 
other values of p, the results are qualitatively and quantitatively the same as that we 
get. 

We thank Prof. Hong Zhao for valuable discussions. This work was partly 
supported by the National Natural Science Foundation of China under Grant No. 
10305005 and the Natural Science Foundation of Gansu Province. 
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Figure 3: Log-log plot of the coefficient F as a function of K. The linear property 
of the plot predicts the power law, F = K~ a . Simulations give the parameter value 
a « 0.98. 
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